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Abstract 

The sum-rate capacity of the parallel Gaussian interference channel is shown to be achieved by 
independent transmission across sub-channels and treating interference as noise in each sub-channel if 
the channel coefficients and power constraints satisfy a certain condition. The condition requires the 
interference to be weak, a situation commonly encountered in, e.g., digital subscriber line transmission. 
The optimal power allocation is characterized by using the concavity of sum-rate capacity as a function 
of the power constraints. 

I. Introduction 

Parallel Gaussian interference channels (PGICs) model the situation in which several transceiver pairs 
communicate through a number of independent sub-channels, with each sub-channel being a Gaussian 
interference channel. Fig. 1 illustrates a two-user PGIC where a pair of users, each subject to a total 
power constraint, has access to a set of m Gaussian interference channels (GIC). Existing systems that 
can be accurately modelled as PGICs include both wired systems such as digital subscriber lines (DSL) 
and wireless systems employing orthogonal frequency division multiple access (OFDMA). Both of these 
systems have been and will be major players in broadband systems. 

While there have been extensions of information theory for the classical single-channel GIC to the 
PGIC (e.g, [1]), most existing research, especially for DSL systems, often relies on the following two 
assumptions [2]-[4]: 
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• Transmissions in sub-channels are independent of each other. 

• Each receiver treats interference as noise. 

These assumptions greatly simpUfy an otherwise intractable problem. The independent transmission 
assumption ensures that the total sum rate is expressed as a sum of all sub-channels' sum rates. The 
assumption of using single-user detection permits a simple closed-form expression for the rate pair of each 
sub-channel ^ Our main goal in this paper is to provide a sound theoretical basis for such assumptions, i.e., 
to understand under what conditions such a transceiver structure leads to optimal throughput performance. 
It is certainly not obvious that this structure could ever be optimal, but we show that it is optimal for 
systems with weak interference, a situation encountered in many deployed systems such as DSL. 

Our approach in characterizing the sum-rate capacity of a two-user PGIC leverages recent breakthroughs 
in determining the sum-rate capacity of the GIC under noisy interference [6]-[8]. We determine conditions 
on the channel gains and power constraints such that there is no loss in terms of sum rate when we impose 
the above two assumptions. This is accomplished in two steps. First, under the independent transmission 
assumption, we find conditions such that the maximum sum rate can be achieved by treating interference 
as noise in each sub-channel. The key to establishing these conditions is the concavity of sum-rate capacity 
in power constraints for a GIC (cf. Lemma 2). Second, we show that with the same power constraints 
and channel gains obtained in the first step, independent transmission and single-user detection in each 
sub-channel achieves the sum-rate capacity of the PGIC. The proof utilizes a genie-aided approach that 
generalizes that of [6]. 

This paper is organized as follows. In Section n we introduce the system model and review recent 
results. In Section m, we consider the maximum sum rate of a special transmission scheme, i.e., 
independent transmission and single-user detection for each sub-channel. We obtain conditions on the 

'Even under this simplified assumption, finding the optimal power allocation is an NP hard problem [5] 
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power constraints and channel coefficients under which the above strategy maximizes the total sum rate. 
We prove in Section IV that the maximum sum rate we obtain is the sum-rate capacity. Numerical results 
are given in Section V. Section VI concludes the paper. 



II. System Model and Preliminaries 
The received signals of the ith sub-channel i = 1, • • • , m are defined as 

Yii = \/ciXii + ^/aiX2i + Zii, 

where < < di,0 < 6^ < q; Z\i and Zii are unit variance Gaussian noise, the total block power 
constraints are P and Q for users 1 and 2 respectively: 



(1) 



and 



E 

i=l 



E 

1=1 



<p. 



where n is the block length, and X\ij and Xnj, j = 1, . . . , n, are the user/channel input sequences for 
the zth sub-channel. We remark that this model is a special case of the multiple-input multiple-output 
(MIMO) GIC [9]. We denote the sum-rate capacity of the ith sub-channel as Ci {Pi,Qi), where Pi and 
Qi are the respective powers allocated to the two users in this sub-channel. 

To find the sum-rate capacity of the PGIC, we need to solve three problems: the first problem is 
whether the sub-channels can be treated separately like the parallel Gaussian multiple-access channel 
[10] and parallel Gaussian broadcast channel [11]-[13], i.e., whether the sum-rate capacity of the PGIC 
is in the form of YliLi^ii^i^Qi)- Such a strategy is suboptimal for PGICs in general [14], [15]. The 
second problem is the optimal distribution of the input signals. It has been shown respectively in [16]-[18] 
and [6]-[8] that Gaussian inputs are sum-rate optimal for a single-channel GIC under strong or noisy 
interference. However, whether this is still the case for PGICs is not known. The third problem is to find 
the optimal power allocation among sub-channels. Existing works on this problem treat the sub-channels 
separately, they use Gaussian inputs, and they use single-user detection at the receivers [2]-[4]. 

Before proceeding, we introduce the following notation. 

• Bold fonts X and X denote vectors and matrices respectively. 

• I denotes the identity matrix and denotes the zero matrix. 
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• |X|, X-^, X~^, denote the respective determinant, transpose and inverse of the matrix X. 

• a;" = [xi,X2, . . . , Xn] ^ is a long vector which consists of the vectors Xi,i = 1, . . . , n. 

• ar ~ (0, means that the random vector x is Gaussian distributed with zero mean and covariance 
matrix T,. 

• E{-) denotes expectation; Var(-) denotes variance; Cov(-) denotes covariance matrix; /(•; •) denotes 
mutual information; h{-) denotes differential entropy with the logarithm base e and log(-) = loge(-). 



A. Noisy-interference sum-rate capacity 

The noisy-interference sum-rate capacity for single-channel GICs [6]-[8] is summarized as follows. 
Lemma 1: The sum-rate capacity of the zth sub-channel with ai < Ci,bi < di and power allocation p, 
q is 



C,{p,q) = \\og{l + -^ , 

2 \ I -\- aiq J 2 



(2) 



1 + hp, 

provided (p, q) ^ Ai'. 

f 

/a^(l + hip) + \/Mi(l + aiq) < \fc^i 

p>0, q>0 

In the case of a symmetric GIC, i.e., = bi,Ci = di and p = q, the noisy interference condition reduces 
to 



A = < 



> . 



(3) 



ai , 1 



< T' P = <1< 



In the case of a ZIC where Oj = 0, the noisy interference condition reduces to 

bi<l, p>0, q>0. 



(4) 



(5) 



The main difficulty in maximizing J2iLi Ci {Pi, Qi) is that Cj {Pi,Qi) is generally unknown if (Pj, Qj) ^ 
Ai- To solve this problem we use the following results. 
B. Concavity of sum-rate capacity 

The key to our study of the PGIC is the concavity of the sum-rate capacity as a function of the power 
constraint. We establish a slightly more general result by using a modified frequency division multiplexing 
(FDM) argument [19]. 

Lemma 2: Let C^(p, q) denote the weighted sum rate capacity of a GIC with powers p and q: 

C^ip, q) = ^ max {Ri + 1J.R2}, 

Ki,H2 achievable 
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where /x > is a constant. Then C^ip^q) is concave in the powers {p.,q), i.e., for any < A < 1 we 
have 

C^^iv, q) > XC^ip', q') + (1 - X)C^{p", q"), (6) 
where p', p", q', and q" are chosen to satisfy 

Ap' + (1-Ay'=p, Xq' + {1 - X)q" = q. (7) 

Proof: Consider a potentially suboptimal strategy that divides the total frequency band into two sub- 
bands: one with a fraction A and the other with a fraction 1 — A of the total bandwidth. Powers are 
allocated into these two sub-bands as {Xp',\q') and ((1 — A)p", (1 — X)q"), where p',q',p",q" are 
such that (7) is satisfied. The information transmitted in these two sub-bands is independent and the 
decoding is also independent. Then the maximum weighted sum rate for the first sub-band is reduced by 
a factor A and becomes XC^{p' ,q'). Similarly, the maximum weighted sum rate at the second sub-band 
is (1 — X)Cf^{p", q"). Therefore, the right-hand side of (6) is an achievable weighted sum rate. ■ 
Lemma 2 provides a fundamental result for weighted sum-rate capacities. It applies not only to two-user 
GICs but also to many-user GICs, Gaussian multiaccess channels, and Gaussian broadcast channels. 

C. Subgradient and subdijferential 

To apply the concavity of sum-rate capacity, we need to use several properties of subgradients and 
subdifferentials (see [20]). 

Definition 1: If / : 7?." — 7?. is a real- valued concave function defined on a convex set S C 7?.", a 
vector y is a subgradient at point Xq if 

f{x)-f{xo)<y^{x-xo), y xeS. (8) 

Definition 2: For the concave function / defined in Definition 1, the collection df (xq) of all subgra- 
dients at point Xq is the subdijferential at this point. 

If the function / is differentiable at xq, then the subgradient and subdifferential both coincide with 
the gradient. We introduce a lemma related to subdifferentials which we use to prove our main result. 

Lemma 3: Let fi{x), i = l,...,m, be finite, concave, real-valued functions on iS C TZ^ and let 
X* e S, i = 1, - ■ ■ ,171. Jf there is a vector y such that y G dfi{x*), i = 1, . . .m, and Yl^i = ^^^^ 
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X* = [xl^, ■ ■ ■ is a solution for the following optimization problem: 

m 



m 

max 

i=l 
m 



subject to ^ Xi =u, Xi e S. (9) 

i=l 

Proof: Since y G dfi{x*), i = 1, . . . m, we have 

fi{x)<fi{x*)+y^{x-xl), Vxg5. (10) 
Let Xi,i = 1, . . . ,m, he any vectors satisfying Xi e S and J2iLi = then using (10) we have 

fi{xi)<Mx*)+y^ixi-x*) (11) 

Therefore 

m m / m m \ 

i=l 1=1 \i=l i=l / 

m 

i=l 

where the last equahty is from Yl'iLi = Z^I^i = ■ 
In the Appendix, we compute the subdifferential dCi{p,q) when {p,q) G A. We are also interested 
in the set of pairs 

Bi= (j dCi{p,q). (13) 
The mapping from Ai to Bi is illustrated in Fig. 2. As seen from (3), Ai is a triangle region with the 



comer points 



and 



0(0,0), 

-5 0, = = (0, qs), 

V airJbidi 



T[ ^^ ,0) ^(pt,0). 



The corresponding points in Bi are respectively 

°' (!■!)■ 



5' 



Ci bidiqs d 



2(1 + Giqs) 2(1 + dig,) 2(1 + d^g,) 
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and 



r 



aiCipt 



2(1 + apt) ' 2(1 + bipt) 2(1 + cpt) 




Fig. 2. The mapping from Ai to Bi. 



Let ^j^^^ be the inner points of At and the Une segment ST, and let B^^^^ be the inner points of the 
closed area defined by O'S'T' and the curve S'T'. As shown in the Appendix, A^^^ maps to B^^^ and 
this mapping is one-to-one. Let A^^ be the line segment OT and B^"^^ be the curve O'T' and the points 

(2) (2) 

above it (labeled as region II in Fig. 2). A] maps to B^ and this mapping is one-to-many. Specifically, 
let {Pi, 0) be a point on OT. The partial derivatives of Ci{p, q) with respect to p and q at this point are 
denoted as Kp and Kq, respectively, where Kp is a two-sided partial derivative and Kq is a one-sided 
partial derivative. Point {Kp, Kq) is on the curve Cyf' of Bf\ The subdifferential of Ci{p,q) at point 

f2) 

{Pi, 0) is a ray in defined as fcp = Kp, kq > Kq. 

Similarly to the above, let be the line segment OS and B^^^ be the curve O'S' and the points 

(3) (3) 

to the right (labeled as region III in Fig. 2). A\ maps to B] and this mapping is also one-to-many. 
Let A^f* be the origin and let B^"^ be the collection of points {kp,kq) satisfying kp > Kp and kq > Kq 
(labeled as region IV in Fig. 2), where Kp and Kq are the two one-sided partial derivatives at the origin. 
A\^^ maps to bI^\ 
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D. Concave-like property of conditional entropy 

The following Lemma is proved in [9] based on the fact that a Gaussian distribution maximizes 
conditional entropy under a covariance matrix constraint [21]. 

T 



Lemma 4: [9, Lemma 2] Let arf 



T T 
1 , . . . , Xj 



, i = 1, . . . , /c, be A; long random vectors each of 



which consists of n vectors. Suppose the Xj^j, i = 1, • • • , all have the same length Lj, j = 1, • • • , n. 
Let = [yf , . . . where has length Lj, be a long Gaussian random vector with covariance 

matrix 

k 

Cov(y") = ^A,Cov(xr), (14) 

where Y^i=\ = > 0. Let 5 be a subset of {1,2,..., n} and T be a subset of <S's complement. 
Then we have 

k 

^ Xih (Xi^s \Xi,T )<h{ys Ivt ) • (15) 

i=l 

When Xk, = 1, • • • , n are all Gaussian distributed, Lemma 4 shows that h {xs \x^ ) is concave over 
the covariance matrices. 

III. A LOWER BOUND FOR THE SUM-RATE CAPACITY 

If the sum-rate capacity of a PGIC can be achieved by (1) transmitting independent symbol streams 
in each sub-chaimel and (2) treating interference as noise in each sub-chaimel, we say this PGIC has 
noisy interference. Before proceeding to the main theorem of noisy-interference sum-rate capacity, we 
first consider the following optimization problem: 



max ^Ci{Pi,Qi 



i=l 

ni in 

subject to ^ = P, ^Qi = Q 

i=l i=l 

Pi>0, Qi>0, i = l,...,m. (16) 

Problem (16) is to find the maximum of the sum of the sum-rate capacities of individual sub-chaimels and 
the corresponding power allocation. In general, the optimal solution of (16) is not the sum-rate capacity 
of the PGIC, since it presumes that the signals transmitted in each sub-chaimel are independent and no 
joint decoding across sub-chaimels is allowed. However, solving problem (16) is important to derive the 
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sum-rate capacity of a PGIC. We are interested in the case where the optimal power allocations P*, Q* 
satisfy the following noisy interference conditions 

V^^(l + 6iPj*) + A/Mi(l + OiQ*) < a/c^, i = l,...,m. (17) 

In such a case, it turns out that the sum of the sum-rate capacities in (16) is maximized when each 
sub-channel experiences noisy interference. 

For the rest of this section, we first consider the general PGIC and derive the optimal solution of 
problem (16) based on Lemmas 2 and 3. We further find conditions on the total power P and Q such 
that the optimal solution of (16) satisfies (17). Then we focus on symmetric PGlCs and provide some 
insights on this solution. 



U 



(18) 



A. General parallel Gaussian interference channel 

Theorem 1: For a PGIC defined in (1), if y^OjCj -|- \/hjdi < y/cidl and the power constraint (P, Q) is 
in the following set 

[k;,klf edCiiP*,Q*) i = l,...,m. 
then the optimal solution of (16) satisfies (17). 

Proof: The proof is straightforward from Lemma 3. For any [k*, k*]-^ G HI^i there exist P*, Q* 
such that [k;,k;f e dCiiP*,Q*),i = l,...,m. Thus if P = Y.T=iPt^Q = YZiQl then from 
Lemma 3, P*,Q* are optimal for the optimization problem (16). Since (P*,Q|) G Ai, then from 
Lemma 1, {P*,Q*) satisfies (17). ■ 
Theorem 1 provides conditions on the power and channel coefficients such that treating interference as 
noise (or single-user detection) maximizes the sum rate of a PGIC under the assumption of independent 
transmission among sub-channels. The conditions of Theorem 1 ensures that the power constraints P and 
Q are associated with a subgradient [k*, k*] shared by Ci{P*,Q*) for alH = 1, • • • , m. Therefore, at 
the points of the optimal power allocations (P*, Q*), all the functions Ci{pi, qi) have parallel supporting 
hyperplanes. We will discuss this in more details in Remark 1 below. 

In general, the closed-form expression (18) of the power region for P and Q is very complex. However, 
for some special cases like symmetric PGICs, we can obtain simpler closed-form solutions. 
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B. Symmetric parallel Gaussian interference channel 

In this section, we consider PGICs with symmetric parameters, namely ai = bi,Ci = di and P = Q. 
Without loss of generality we assume ci > C2 • • • > c^. Define 



w = maxjtUj}, 



2 ' 



and let r to be an index between 1 and m such that 



Cr+l < W < Cr, 

where we let Cm+i = for convention. Then we have the following theorem. 

Theorem 2: For a symmetric PGIC, if — <\, i = 1, . . . ,m, and 

Ci 4 



(19) 
(20) 

(21) 



< P < P, 



(22) 



where 



E 



i + — ^(«» + Cj) - (2aj + Ci) 



2ai{ai + Ci) 

t — J. 

then the optimal solution of (16) satisfies (17). Furthermore, only the first r sub-channels are active. 
Proof: By symmetry, we simplify the proof by considering the following optimization problem: 



(23) 



max ^Q(Pi,Pi) 

i=l 

rn 

subject to ^ Pi = P 



(24) 



i=l 



Pj > 0, i = 1, . . . ,m. 

That is, we require that the power allocated to both users be Pj for the ith sub-channel. Obviously the 
maximum of (24) is no greater than the maximum of (16) because of the extra constraint Pj = Qi. To 
prove Theorem 2, it suffices to show that under the condition < P < P: 1) the optimal P* for (24) 
satisfy the noisy interference condition; 2) the optimization problems (16) and (24) are equivalent. 
Let Pi = Qi, we obtain from Lemma 1 



P 



0<p< 



2ai 



(25) 



The subdifferential is computed in the Appendix and is given by (see (89)) 

{k\ci<k<c}, Pi=0, 



dCi{Pi] 



k 



k 



{l + aiPi){l + aiPi + CiPi_ 



, PieX, Pi 7^0, 



(26) 
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where c = max{ci}. Therefore 

I3'i= U dCi{Pi) = {k\wi<k<c}, (27) 

and 

m 

f]Bl = {k\w <k<c}. (28) 

i=l 

For any k* G 07=^1 ^'v equation (26) determines a one-to-one mapping from k* to P* G .4^, namely 

0, k* >Ci 

^*(^*) = ,/c| + ^(a, + c.) - (2a. + c.) (29) 



2ai{ai + Ci 
So consider the region (18) which is here 



Wi < k* < Ci 



U 

fe*e[t&,c] 



P = J2Pt{k*)\- (30) 
i=i J 

From (29), X^^^ is decreasing in A;*, therefore 

p > p{k* = c) = 0, (31) 
P < P{k* = w) = ^P*{k* = w) 



i=l 



'cf + ^^{ai + Ci)- {2ai + 
^ 2ai{ai + Ci) 

where the first equahty of (32) is from (21). Since P(A;*) = Y^iLxPti^*) is continuous over k* , for 
any P G [0, P] there exists a k* , and the corresponding P?',i = 1, . . . ,rn, that solve the optimization 
problem (24). 

We complete the proof by showing that the optimal P* for (24) is also optimal for (16) for a symmetric 
PGIC. Assume that for a given P and the optimal P* of (24), the corresponding subgradient (which is 
identical for all i) is k* . Then by symmetry, the subderivative of Ci{Pi, Qi) in (16) is by choosing 

Pi = Qi = P*- Therefore the subderivatives are identical for all the Ci{Pi, Qi) at Pi = Qi = P*. From 
Lemma 3, Pi = Qi = P* is an optimal choice for (16). Since P* satisfies the noisy-interference condition 
in (25), Pi = Qi = P* also satisfies the noisy-interference condition in (3). ■ 

Remark 1: From the proof of Theorem 2, all the Ci{Pi, Qi) have parallel supporting hyperplanes at 
the optimal point Pi = Qi = P*. This gives rise to a geometric interpretation, as illustrated in Fig. 3. 
For clarity, we use the simplified optimization problem (24). Ci{p) is the sum-rate capacity for the ith 
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Fig. 3. An illustration of sum-rate capacity achieving power allocation for a symmetric parallel Gaussian interference channel. 



sub-channel, points Ai and A2 correspond to the power allocations Pi and P2, respectively, and the two 
supporting hyperplanes pass through Ai and A2. The power allocation satisfies P = Pi+P2,Pi G A'^,i = 
1,2, and k = ^'^gjf^ \p=Pi = ^^djf^ \p=P2 (we assume the subgradient is equal to the gradient in this 
case, and hence the supporting hyperplane is the tangent hyperplane), and the corresponding sum rate is 
Cli + C*2- Consider the power allocation Pi— 6 and P2 + S and the corresponding sum-rate capacities for 
the two sub-channels C*i — ARi and C*2 + Ai?2> respectively. By concavity, Ai?i > k6 and AR2 < k6. 
Therefore the new sum-rate is {C*^ - ARi) + {C*2 + AR2) < {C*^ - k6) + {C*2 + k6) = C*^ + C*2. 

Remark 2: When min | | < P < P, there exist power allocations such that some sub- 

channels do not have noisy interference. As such the sum-rate capacities of those sub-channels are 
unknown. Surprisingly, in this case we do not need to derive upper bounds for those unknown sum-rate 
capacities. Instead, the concavity of the sum-rate capacity (as a function of the power) and the existing 
noisy-interference sum-rate capacity results ensure the validity of Theorem 2. 

Remark 3: The parallel supporting hyperplanes condition for the optimal power allocation is applicable 
to a broad class of parallel channels in which 1) transmissions across subchannels are independent, 2) 
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the capacity of each subchannel is concave in its power constraint. For example, this condition apphes to 
parallel multi-access and broadcast channels. In particular, applying the condition to single user parallel 
Gaussian channels, it is easy to verify that the parallel supporting hyperplanes condition reduces to the 

classic waterfilling interpretation. 

Remark 4: Intuitively, since each sub-channel is a symmetric Gaussian IC with noisy interference, the 
power allocated to the two users in each sub-channel ought to be identical. While Theorem 2 does not 
explicitly address the power allocation scheme, we see from the proof of the theorem that this is indeed 
the case. 

Remark 5: If (22) is satisfied, the optimal power allocation P* is unique, and there exists a A;* G [it), c] 
such that P* and k* satisfy (29). To see this, observe that in the proof of Theorem 2, YhLi P* {k*) is 
continuous and monotonically decreasing in k* when k* G [w,c], and P varies from to P when k* 
varies from c to w. Thus, if < P < P, there exists a corresponding unique k* in [w, c] that solves 
problem (24). Since the mapping from k* to P* in (29) is a one-to-one mapping, P* is also unique. 

Remark 6: As shown in (29), whether a sub-channel is active or not depends only on the direct channel 
gain Ci. The amount of power allocated to a sub-channel depends on both the direct channel gain Ci and the 
interference channel gain a^. When the total power constraint P increases from to P, the corresponding 
k* decreases from c to w. As such, from (29) the sub-channels with larger Cj become active earher than 
those with smaller q. 

IV. Noisy interference sum-rate capacity 

The following theorem gives the noisy-interference sum-rate capacity of a PGIC. 
Theorem 3: For the PGIC defined in (1), if 



for all i = I,-- - ,m, and the power constraint pair {P,Q) is in the set (18), the sum-rate capacity is 
the maximum of problem (16), and the sum-rate capacity is achieved by independent transmission across 
sub-channels and treating interference as noise for each sub-channel. 
The following theorem is a special case of Theorem 1 for symmetric PGICs. 



satisfies (22), then the sum-rate capacity is the maximum of problem (24) and is achieved by independent 
transmission across sub-channels and treating interference as noise in each sub-channel. 




(33) 



Theorem 4: For a symmetric PGIC, if — < - for all ? = 1, • • • 



m, and the power constraint P 
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For the PGIC, there may exist some sub-channels with one-sided interference or no interference, i.e., 
bj = 0, or Ofe = or Or- = by = for some integers j, k and r between 1 and m. We prove Theorem 3 
for all such cases. 

Proof: Let i,j,k,r be integers and 1 < i < mi, mi + 1 < j < m2, m2 + I < k < m^ and 
777,3 + 1 < r < m throughout this proof. We denote i,j,k,r as index sets and i = {!,••• ,mi}, 
j = {mi + 1, ■ • • , 771,2}, k = {m2 + !,••• , m^}, and r = {m^ + !,••• , m}. Without loss of generality we 
can assume that sub-channels with index i are all two-sided GICs with Ui ^ 0, bi ^ ; the sub-channels 
with index j are all GICs with aj 7^ 0, bj = 0; the sub-channels with index k are all GICs with Uk = 0, 
bk 7^ 0; and the sub-channels with index r are all GICs with Ur = br = 0. Let A = diag(Y^, • • • , ^/a^), 
B = diag(v^, • • • , Vb^), C = diag(y^, • • • , v^)» D = diag(v^, • • • , Vd^)- Then we can rewrite 
(1) in the following form 

yi = Cxi + Ax2 +zi, 

(34) 

y2 = Bxi + 'Dx2+Z2, 

where zi ^ Af (0, 1) and Z2^ M (0, 1). 
We further define 

Aj = diag(^/ar, • • • , v'a;;^) = diag(^), 
Aj = diag(^a^j+i, • • • , ^/a:^) = diag(^/^), 
Afc = diag(^am2+i, • • • , ^/a^) = diag(^/a^), 
A^ = diag(Va„3+i, • • • , ^Ja^) = diag(y^), 

and similarly for B, C and D. Denote the transmitted vector of user 1 as Xi = [Xi,- ■ ■ , X^], where each 
entry is the transmitted signal at the corresponding sub-channel. Similarly, we let Xii = [Xi, ■ ■ ■ , X^J = 

[Xli], Xlj = [Xrrii+l, ■ ■ ■ , Xrn^ = \Xij\, Xik = [Xm2+1, • ' ' , Xm.^] = [Xik] andXlr = [Xjus+l, ■ ■ ■ , Xm] = 

[Xir]. The input vectors for the second user are similarly defined. 
Since the power constraint [P, Q]'^ is in the set (18), so there exists a subgradient [k*, k*]^ G H/^i 

rp 

and the corresponding [-Pf , ^ -^l such that 

[k;Xf &dCi{Pt,Qr), l = l,---,m. (35) 
From Theorem 1, the P^*, optimize problem (16). 
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Assuming the channel is used n times, the transmitted vector sequences are denoted as a;" = [xj^ , ■ ■ 

T 



J' -1 T 



and X2 = [x2i , • • • ) ^2n\ which satisfy the average power constraints 

n 

Y,^x[E{xnxl)] <nP, 

1=1 

n 

J2^r[E{x2ixl)] <nQ. 



1=1 

Define zero-mean Gaussian vectors x* = ^X*i , ■ ,X 
matrices 



Im 



andxo = 



■^21 1 " ) 



2m 



1 " 

Cov (xi) = - ^ Cov {xii) , 
^ 1=1 

1 " 

Cov (^2) = - XI ^^'^ *^^20 • 



Obviously, and X2 satisfy the power constraints. We define 

P^=Var(xr,), 
Qi = Var (X*2i 



with the covariance 

(36) 
(37) 

(38) 
(39) 



Vectors yl and are defined by (34) with xi and X2 being replaced by xl and X2 respectively. Similar 
to xi, xl is also partitioned as x^j, x^j, x^i^ and xl^.. 

Define Gaussian random vectors nii,nij,n2i and n2k independent of Xi and X2, and let 







( 


h 








Plify li 





\ 


Zij_ 









h 
















0, 








Ifc 










nii 






Plicr 1, 








-I 
















Plj(Tlj_ 








. 


/ 


Z2i 




/ 


h 








P2iO'2i 





\ 


Z2j_ 






















Z2k 




0, 








\ 





P2kf^2k 




n2i 






P2iO'2 







4 







n2k_ 




V 





P2k<^2k 





^Ik . 


1 



(40) 



(41) 
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where pi and ai are diagonal matrices with the diagonal entries being pi and (jj, i E i, respectively. 
Furthermore, we let 



a 



li 



4 



1 



^ (a,Q* + 1)' - ^ m + 1)2 + 1 



4/) 

— {aiQ* + 1)^ 



«^ (6,i^* + 1)2 - ^ (a,g* + 1)2 + 1 



P2i 



Pli 



2 



P2fc 



± 



+ 1)2 _ ^ (a-Q* + 1)2 + 1 



Cfc' 



(42) 



(43) 

(44) 
(45) 

(46) 
(47) 
(48) 

(49) 



We emphasize that the and Qi in (42)-(49) are the optimal powers for the problem (16) and can 
be considered as constants in what follows. It has been shown in [6, equations (51),(52)] that (42)-(45) 
are feasible (i.e., there exist at least one choice of {crfj, crfji Pii' P2i} such that the covariance matrices 
are symmetric and semi-positive definite, and thus the defined Gaussian random vectors exist) for the 
definition in (40) and (41) if and only if [P*,Q*f G A- Obviously, (46)-(49) are feasible for the 
definitions in (40) and (41) if and only if aj < dj and bk < Ck- Moreover, (42)-(45), (46) and (48) satisfy 



Vcipiiaii = 1 + aiQi, 
Vdip2icr2l = l + biPi*, 



(50) 
(51) 



for alH = 1, • • • ,m. 
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Let e > and e — >^ as n — oo. From Fano's inequality, any achievable rate Ri and R2 for the PGIC 
must satisfy 

n{Ri + R2) — ne 

< I (xl-y^xl + nl_, x\ + <,) + / (A-.vl.Ai + + «2fe) 

= h {y^,xl + nl,xl-+nl^ - h +n?,,x?^ + | a:?) + h {y^,xl^ + nl,x^,+n^,) 

-h {y^,xl_ + nl,x^, + n^fc I ) . (52) 

In (52), we provide side information x"j and x^j to receiver one, and X2i +«2i ^2k +'*2ifc 
receiver two, respectively. For the first mi sub-channels which have two-sided interference, both receivers 
have side information. For the sub-channels which have one-sided interference, only the receivers suffering 
from interference have the corresponding side information. For the sub-channels without interference, no 
side information is given. 
For the first term of (52), we have 

7 f .,n n n n n , n n , n \ 

— " I i/liJ i/lj) i/lfe) 2/lr ) -"^li ~r "li) -Cij + 11^] 

< h {y'lk,<: + «h) + h (yl \ x\ + n?, ) + {fi^_,x\ + + (y?,) 

< h (y\-<^\ + nn) + | + ) + nh (yl^_, x^,^ + ) + nh {y\^}j 

< h (y^, , xl +nl^+n^h (y^^ | X^, + iV^ ) + n ^ , X^,- + iVi,) + n ^ /i (y^;) ,(53) 

i j r 

where the first inequality follows by the chain rule and the fact that conditioning does not increase 
entropy, and the second inequality is from Lemma 4. 
For the fourth term of (52), we have 

-h[y2, + n2i, X2k + n2k \ ^2 j 

= —h (^y'ii, V2j ) V2ki V2rJ ^'ii + '^24 ) ^2k + '^2fe ^2i' ^2j ) ^2fc) ^2r ) 
= — /l ^BjX^j -I- «2i) ^2j ) Bfe^lfe + ^2fe) ^2r ) '*2i ) '*2fe) 

= -/i (b,x?, + BfeX^fc + ^r^,, nl,n^2^ - h {z^^ - h (z^,) 

= -h {b,xI + z% B,x^, + z\ I n|, n^, ) - /I (n|) - h (n^,) - (z^^) - h [z^^ 
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= -h {BiXl + zl, BkX^, + z^, I nl, n^, ) - n ^ /i (iV2,) - n /i {N2k) -nj^h {Z20) 

i k j 

—n'^h{Z2r) ■ 

r 

where the third equaUty holds since and z^.^. are independent of all other variables. 
Combine the first terms of (53) and (54), we have 

h {y% x\ + <,) - h {^ixl + zl, Bkx'ik + z\ I nl, 



(a) 



h (xl_ + nl_, x\ + C, iz^,) - h (xl_ + B-^wl, x\ + B^i«,5,) + n log ■ 

'Cfcl 



(6) 1 
= niog 



(c) 



-n^h (xl^ + -^W2k] +n^log^-n^log V^-n^logA/6fe 

i A; j 

k 

where in (a) we let W2i and tW2fc be independent Gaussian vectors and 

Cov {w2i) = Cov {z2i \n2i) =li_- diag(p2j) 
Gov (w2fe) = Gov {z2k \n2k) =lk- diag(pife). 



(54) 



(55) 



(56) 



The stacked vectors and W2fc each have independent and identical distribution (i.i.d) entries. Equality 
(b) holds because of (45) and (49) which imply 

Gov (nii) = Gov ^B~^t(;2i) , 
Gov (C^'^ife) 

and 



Gov (B^^t/)2fe) , 



h (xl_ + x^, + C^'z^,) - h (xl + BT'wl, x\ + B-^wl^ = 0, 
regardless of the distribution of a;"j and a;"^. 
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Equality (c) holds also because of (45) and (49), which imply 

h (li*, + iVn) - h + -^W2^^ = 0, 
h f + -^Z,k] - h (xlk + -^W2^ = 0, 



Ir 



regardless of the distributions of X\i and Xi^. 
Combining (53) and (54) and using (55), we have 

h {yl,xl_ + nl^ , xl^_ + 71%^ - h {y^,xl + n% x\ + n\ \ 

i A; i »■ 

-n ^ (iVaO - n ^ (iV2fe) - n ^ /i (Z2,) - n ^ (Z2.) . (57) 

Similarly, because of (44) and (47) we have 

h {f2.x\^n%xl^-rnl^ - h {y'l,xl+nl,xl_ + nl_ \ x?) 

i j i r 

k i j 

-nY,h (Nu) -nY,h{Nij)-nY,h (Z^k) -n^h (Zi,) . (58) 

i j k r 

Substituting (57) and (58) into (52), we have 
Ri + R2 — € 

<Y,[h (X*, + Nu) + h (y* I X;, + Nu)-h (Nu) - h [v^iX*^ + Zu \ Nu ) 

i 

h (^2* + N2i) + h (9* I X*, + N2i)-h {N2i) - h {VhiX*u + Z2i I N2^ 

+ E [^^ i^yv% + ^li) - ^ (^li) - ^ (V^^2,- + Z^j I iVi,- ) + h + ^2,) - h {Z2j) 

3 

+ (^2fc, ^2fc + N2k) - h {N2k) - h [^kXlk + Z2k \N2k)+h {yTkXlk + Zik) - h {Z^k) 

k 

+ E (n;) - h {Zu)\ + ^ [/I (^2;) - h (^2.)] 

r r 

i j k r 
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1=1 

where 

1 



(59) 



fii-^ij Qi) ~ 






+^log 


fj{Pj,Qj) 




fk{Pk,Qk) 





li 



1 + aiQi - \aii 1 + QiQi 



{l + biPi)Qi ( 1 V^p2 



+ 1 + 



1 + aiQi 



il + ajQj)Pj f 1 ^pij 



1 + OjQj — PI ■ Vc^ij 1 + c-jQj 



+ 1 + 



+ 1 + 



(1 + bkPk) Qk (_}__ \/4p2fc y _|_ ^ _|_ 4Qfc 

1 + 6fcPfc - pifc Vc^2fc l + bkPkJ l + bkl 



P2k 

fr{Pr, Qr) = \ log (1 + Pr) + \ log (1 + Qr) ■ 



1 + hPi 
1 + OjQj 

1 + ^'ifcf'fe 



(60) 



+ - log {I + djQj),m 



+ - log (1 + CfePfc) , (62) 



(63) 



Next we will show that the fi{Pi,Qi), I = 1, • ■ ■ are all concave and non-decreasing functions of 
{Pi, Qi). From (59) and the fact that 

h (Xli + iVi,) - h {^fbiXl, + Z2i I iV2i ) = - log v^, 

/i (xli + A^2i) - ^ ( Vai^2*i + Zu \Nu)=- log ^H^, 

we have 

fi{Pi, Qi) = h f Fi* XI, + Nu)- h (Ni,) + h (%\ XI, + N2^\-h {N2i) - log (64) 



Define Gaussian variables X^- ^, X^^ ^ and ^ independent of Ni, and Yj* ^ = y^^ii t+\/^^2i t+-^i' 
t = 1, - ■ ■ ,s where s is an integer. Let Var (^X^- = P^j, Var = Qj j and 



5^ AtVar = Pi = Var (Xi*,) , 

t=i 

s 

At Var (X^V,) = = Var (x^*,) , 



(65) 
(66) 



where {Aj} is a non-negative sequence with X^^^^ -^t = 1- Then we have 



^A^Cov 



t=i 



Gov 



(67) 



From Lemma 4 we have 



h (n* I Xli +Nu)>Y\th (nVt I ^lU + ) 



Z=l 



(68) 
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Therefore, h ( Y^- 



N2i) 



-^li + -^li ^ is a concave function of (P^ , Qi). For the same reason h (^2i 
is also a concave function of {Pi,Qi). Therefore fi{Pi,Qi) is a concave function of {Pi,Qi). Similar 
steps show that fj, fk, and fr are concave in {Pi, Qi). 

To show that fi{Pi,Qi) is a non-decreasing function of Pi,Qi, we let Xu, X2i, Ui and U2 be four 
independent Gaussian variables and 

^li = + Ui, 

^21 = + ^^2- 



Then 

|xi*, + Arn) 

= h (Vq^h + V^^2*. + Zii \xii + iVii) 

< h (v^Xi*, + V^X2*, + Zu \x;, + ATu, C/i , C/2 ) 

= ^ (-/q-^h + ^/aiX2^ + Zii \Xii + Nii) . 

Therefore, h ^Yj* X*^ + Nu ^ is a non-decreasing function of {Pi,Qi). For the same reason, h (Y2i 
is also a non-decreasing function of {Pi,Qi). Therefore, fi is a non-decreasing function of {Pi,Qi) 
From (50) and (51) we have 



*2^+N2i) 



m ^ m 



=1 



log 1 + 



ciPr 



1 + aiQ* 



+ log 1 + 



i + hp; 



(69) 



1=1 



Next we will show that YZif{Pi*,Ql) > TZifiPhQi) for any Pi,Qi that satisfy TZiPi = P' 

TT=iQi = Q- 

Using (50) and (51), we have 



dp 

dfiip, 



dq 



i = Qt 



p = p' 



dCiip,q) 
dp 

dCiip,q) 
dq 



p = Pt 
i = Qt 



q = Q*l 



(70) 



(71) 



for all Z = 1, • • • ,m. Therefore, // and Ci have the same partial derivatives at point (Pf, Q\). From the 
Appendix, the subgradient of a function is determined by the derivatives, therefore, // and Ci have the 
same subgradient at point {P*,Qf) for each I. From (35), we have 

[k;,k;]edfi{Pt,Qt), l = l,---,m. (72) 
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Therefore, from Lemma 3 we have 

m m 

^fiiPi*,Ql)>Y.fiiPi'Qi)' (73) 
1=1 1=1 

if EI^i^/ = ^and EZiQi = Q- 

If Yl]^iPi ^ P ai^d Y2i=iQi ^ there exist two non-negative sequences {si,-- - ,s„i} and 
{^1) • • • ,tm} such that Ylb=ii-^i + si) = P and Yl^iiQi + = Q- Therefore we have from Lemma 3 

m m 

J2 fi ^PtM) >J2fi {Pi + si, Qi + ti) . (74) 

1=1 1=1 

Since fi{Pi,Qi), I = 1, - ■ ■ , m are all non-decreasing functions of Pi and Qi, we have 

m m 
i=l 1=1 

Combining (74) and (75) we have 

m m 
1=1 i=l 

for any Y^"l^ Pi < P and XlJ^i Qi < Q- Therefore we have from (59) and (69) that 



^ m 

Ri+R2-e<-Y^ 



2 

1=1 



log 1 + . , v.* + log 1 + 



(77) 



The above sum rate is achievable by independent transmission across sub-channels and single-user 
detection in each sub-channel. Therefore (77) is the sum-rate capacity of the PGIC if the power constraints 
satisfy (18). ■ 
Remark 7: The main idea of the proof can be summarized as follows. We first assume an arbitrary 
power allocation (P;, Qi), I = 1, - ■ ■ ,m. Then we show that the sum rate for this power allocation is 
upper bounded by J2iLi fi {Ph Qi)- This upper bound decomposes the sum rate bound into the sum of the 
individual sub-channel's sum-rate capacity upper bounds. By Lemma 3, the maximum of Y^^i fi {PiiQi) 
is YliLi Ci {P^ ,Q*i) which is an achievable sum rate for a special power allocation. To ease the proof, 
the upper bound /; can not be arbitrarily chosen. Compared to the sum-rate capacity Ci for sub-channel 
I, ft has the following properties: 

• /; is concave over the powers; 

• fi is tight at the optimal point (Pf , Q^); 

• /; and Ci have the same subdifferentials at the optimal point {P^,Qi). 

Therefore we choose the noise vectors ni and n2 such that the above conditions are satisfied. Fig. 4 
illustrates such an upper bound. 
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Fig. 4. An illustration of the upper bound suitable for the proof of Theorem 1. 

V. Numerical examples 

Figs. 5 and 6 are examples for a symmetric PGIC with directly link gain ci = di = 1. Fig. 5 shows 
the ratio of maximum noisy interference power constraint P for a two-channel symmetric PGIC with a 
varying from to 0.25, and the sum of the maximum noisy interference power constraints for Si + 5*2 
where Si = ^^^i"" , i = 1,2. Thus, for ai = 02 (i.e., these two sub-channels are identical), the ratio is 
1 and Si = S2 = ^, and when ai and 02 are far apart, then P Si + S2- Thus, for ai = 02 we achieve 
capacity despite the fact that we do not know the capacity if all power is placed in one sub-channel. 
Again, this is where Lemma 2 is useful. 

Fig. 6 shows the maximum noisy interference power constraint P for a two-channel symmetric PGIC 
with ai varying in [0, |] and 02 = |. When ai = j, the first sub-channel no longer has noisy interference, 
therefore the maximum noisy interference power constraint is P = 0. Fig. 6 also shows that P decreases 
with oi. The discontinuity at oi = | is because the second sub-channel becomes the worse channel. 

Fig. 7 shows the regions of Bi and B2 for a PGIC with two sub-channels. In this case Bi C B2- 
In Fig. 8, OiCMND is the noisy-interference power region for this PGIC. Points A,E coincide in 
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Fig. 8 since in Fig. 1 A,E & B^^ [\ and thus the power allocation is Q* = P| = = and 
Pf > 0. Similarly, points C,T coincide in Fig. 8 since C,r G Bf^ {^Bf' and thus = Q* = 
and > 0, P2 > 0- Similar arguments apply to points P, F and P, S*. By considering the relationship 
between the power regions and the respective subdifferentials we can determine the activity of the two 
users in each sub-channel. This activity is summarized in Tab. I, where indicates inactive (zero power) 
and + indicates active (positive power) for the user in the corresponding sub-channel. In the following 
we illustrate the regions in Figs. 7 and 8. 

We first remind the reader that B^^^ corresponds to the case where both users are active; B^^ corre- 
sponds to the case where only user 1 is active; B\ ' corresponds to the case where only user 2 is active; 
and B^^ corresponds to the case where both users are inactive in sub-channel i. Consider the following 
regions. 

• The region O2MN in Tab. 1 denotes regions in both Figs. 7 and 8. In Fig. 7, it is the intersection 
of B^^ and b'^\ Thus, to achieve the sum-rate capacity both users are active in both sub-channels. 
The corresponding power region O2MN is shown in Fig. 8. 

• Region O1EO2F is the intersection of B^^ and B^\ So both users are active only in sub-channel 
1. In this case, both of the power constraints P and Q are small, so that only the better sub- 
channel which produces larger sum-rate capacity than the other is allocated power. Therefore, this 
two-channel PGIC behaves as a GIC. 

• Region OiAE is the intersection of B^^^ and B^\ So user 1 is active in sub-channel 1 and user 2 
is inactive in both sub-channels. The power region 0\AE is shovm in Fig. 8. The overall power for 
user 2 is Q = 0. 

• Similar to the above case, region AETC of Fig. 7 is the intersection of B\ and B2 ■ The optimal 
power allocation makes user 1 active in both sub-channels while user 2 is inactive in both sub- 
channels. In Fig. 8, the overall power for user 2 is also (5 = 0. Actually, regions 0\AE and AETC 
are examples of single-user parallel Gaussian channels whose optimal power allocation is the water- 
filling scheme. In the former case, the power constraint is so small that only the sub-channel with 
larger direct Unk channel gain (sub-channel 1) is allocated power. In the latter case, the power 
constraint increases to a critical level (point A in Fig. 8) so that both sub-channels are allocated 
power. 

• Region ETMO2 is the intersection of and B^2^ . User 1 is active in both sub-channels while 
user 2 is active only in sub-channel 1 which has larger direct channel gain for user 2. In this case. 
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the two-channel PGIC with two-sided interference works like a PGIC with one sub-channel having 
two-sided interference and the other having one-sided interference. 
. Regions OiBF, BFSD and FO2NS are counterparts of regions OiAE, AETC and ETMO2, 

respectively, by swapping the roles of the two users. 
Also plotted in Fig. 8 are the noisy-interference power regions for the individual sub-channels, where 
sub-channel 2 has a larger noisy-interference power region than that of sub-channel 1. In this case, the 
overall noisy-interference power region is larger than that of either of these two sub-channels. 




Fig. 5. Ratio of the total power constraint and tlie sum of tlie individual sub-channel power constraints for different channel 
gains. 



VI. Conclusion 

Based on the concavity of the sum-rate capacity in the power constraints, we have shown that the 
noisy-interference sum-rate capacity of a PGIC can be achieved by independent transmission across sub- 
channels and treating interference as noise in each sub-channel. The optimal power allocations have 
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0.05 0.1 0.15 0.2 



Fig. 6. The maximum noisy interference power constraint P for ai with 02 = |- 



TABLE I 

Power constraints and the activeness of users 





OiA{E) 


A{E){T)C 


ETMO2 


01E02F 














(+,0) 


(+,0) 


(+,+) 


(+,+) 


{P-LQi) 


(0,0) 


(+,0) 


(+,0) 


(0,0) 




02MN 


OiB{F) 


B{F){S)D 


F02NS 












{Pi,Ql) 


(+,+) 


(0,+) 


(0,+) 


(+,+) 




(+,+) 


(0,0) 
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Fig. 7. Bi and S2 for the parallel Gaussian interference channel with a\ = h\ = 0.6, c\ = d\ = 4, 02 = ti2 = 0.24, C2 
d2 = 1.2. 



the property that the sub-channel sum-rate capacity curves have parallel supporting hyperplanes at these 
powers. The methods introduced in this paper can also be used to obtain the optimal power allocation 
and capacity regions of parallel Gaussian multiple access and broadcast channels [22]. 

Appendix: subdifferential of Cj(p, q) 

The subdifferential dCi{p,q) depends on the location of the point {p,q). If the sub-channel is a two- 
sided GIC, from Lemma 1, is defined in (3). We derive the subdifferentials as outlined in (78)-(88) 
below and present the evaluations in (89) below. 

• If {p, q) is an interior point of Ai, Ci{p, q) is differentiable. From the concavity of Ci{p, q), "^Ci{p, q) 
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Fig. 8. Noisy-interference power region for the parallel Gaussian interference channel with ai = 61 = 0.6, ci = di — 4, 02 = 
62 = 0.24, C2 = d2 = 1.2. 



satisfies (8) and tlie subdifferential dCi{p, q) consists of a unique vector [kp, kq\^ = vCi{p, q) where 

dCi{p, q) 



k 



k, 



^ dp 
dCi{p, q) 



dq 



(78) 
(79) 



If ^yaiC^{l + bip) + y/bidi{l + aiq) = ^/cidi,p > 0, g > 0, we can compute only one-sided partial 
derivatives since Ci{p, q) is unknown on the other side. The subdifferential includes a vector [kp, kg^ 



where 



kp = lim 



ko = lim 

where lim = lim and similarly lim = lim 

(5tO (5<0,5-^0 510 (5>0,<5^0 



Ci{p + 6,q) - Ci{p,q) 
510 5 ' 

Ci{p,q + 6) - Ci{p,q) 



(80) 
(81) 
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• If g = 0,p > 0, we have only a one-sided partial derivative in q. From the concavity of Ci{p,q), 
all the [kp,kq]'^ vectors satisfying following conditions are subgradients 



< p < 



dCiip,0) 

dp ' - - bi^/a^ 

,. Ci{p + 6,0) - Ci{p,Q) Vcidi - Ja^Ci - y/b^i 

hm — ^ ^ -, p= , , 

5T0 Oi^aiCi 



ka = lim 



a{p,^)-Ci{p,o) 



(82) 



(83) 



On the other hand, for the same point {p, q) and the corresponding [kp, kqY defined above, by 
choosing k'^ > lim^j^o ^^^^'^-^ ^ Ci{p,0) ^ vector [kp, kq]^ satisfies (8) for all points [p, g]^ G Ai, 
and thus is also a subgradient. Therefore we can replace (83) with 

Ci{p,d)-Ci{p,0) 



kq > lim ■ 



5^0 S 

Similarly, if p = 0,q > 0, the subdifferential is the set of [kp, kq]'^ with 

ai5,q)-aio,q) 



kn > lim 



< g < 



/ djCi 



sio 6 
( da{0,q) 
dq ' 

Ci{p + S,0) -Ci{p,0) ^Jc^i - Ja^i - \/h^i 

lim , q = 

L ^TO 6 aiy/bidi 

If p = g = the subdifferential is the set of [kp, kq]'^ with 



^hdi 



kq > lim . 
For completeness, we summarize dCi {p, q) as follows 



(84) 



(85) 



(86) 



(87) 
(88) 



dCi{p, 



March 3, 2009 



DRAFT 



30 



[kp, kq) 
{kp, kq) 

{kp, kq) 
{kp, kq) 



1 



hi 



2 V 1 + Cjp + ttjg 1 + bip + diq 1 + bip^ 



a,; 



+ 



di 



a,; 



, ip, q) e A 



(1) 



2 \l + Cip + aiq 1 + bip + diq 1 + Uiq^ 

Ci 



2(1 + dp) 
di ttiCiP 



d 

2(1 + bip) 2(1 + dp) - « - 2 

Ci bidiQ c 

< kp < - 



ip, q) e A 



(2) 



2(1 + Giq) 2(1 + dig) 

k = ^» 
« " 2(1 + 

-^<kp<-, -^<kq<- 



(p, q) e A 



(3) 



(p, q) e A 



(4) 



where 



Af^ = {(p, v^(l + bip) + y/bi{l + Giq) < \fddi,p > 0, g > 1 , 



A 



(2) 



A) 



(3) 



(P, 9) 



<p< 



P = 0,0 < g < 



Cy'C^i \/ CLi Ci 



'bidi 



bi-\^^ttiCi 
\J Cidi ttiCi 



-,q = 0}, 
^bdi. 



aiVbidi 



A^^ = {{p,q) \p = q = 0}, 



and 



c = max {cj}, 

i=l,--- ,m 

d = max {di}. 

1=1, ••• ,m 



(89) 



(90) 
(91) 

(92) 
(93) 

(94) 
(95) 



In (89) Ai = Uf=i>^F^- For the subgradient [kp^kq^ , when p = Q ox q = 0, kp or kq varies from 
some constants to infinity. In (89) we upper bound kp and kq with | and | respectively for convenience 
and without loss of generality. The main reason is that we are interested only in IJ™ ^ Bi. To relate the 
mapping of A^p to different regions of Bi, we further define 



b\^^= U dCi{p,q), i = l,---,4. 



(96) 
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If the sub-channel has one-sided interference 6^ = 0, < Oj < 1 or no interference Ui = bi = 0, then 
from Lemma 1 we have Ai = {p > 0,q > 0}. We similarly obtain dCi{p, q) as follows: 



dCiip, q) 



= < 







< 












< 


(fcp, kq) 




{kp, kq) 







^ 2(1 + CiP + Giq) 
k„ = l( ^ . + 



2\l + Cjp + ajg 1 + dj^ l + a^q^ 
^^=2(1 + CP) 

2 2(1 + Qp) - - 2 
c,- , c 

2(l + aig) - P - 2 



di 



2 l + di 



-^<kp<-, ^ < /cq < - L 



p > 0, > 



p > 0, g = 



p = 0, gr > 



p=q = 0. 



(97) 
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